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Abstract 

We use the separation of variables and T-matrix methods to calculate the optical 
properties of homogeneous spheroids with refractive indices from m = 1.3 + 0.0* up 
to 3 + 4i. 

It is found that the extinction cross-sections for highly absorbing spheroids are 
normally 1.5 - 2 times larger than those for spheres of the same volume. The albedo 
of the non-spherical particles rather slightly depends on the particle shape and is 
mainly determined by the imaginary part of the refractive index. Beginning at some 
size, the spheroidal particles do not polarize the transmitted radiation independent 
of their shape. 

We also suggest a new approach for axisymmetric particles which would combine 
the strong aspects of both methods mentioned above and give several values of the 
cross-sections as benchmarks in tabular form. 

Key words: Light scattering; Non-spherical particles; Polarization; Benchmarks 



Accepted for publication in Journal of Quantative Spectroscopy & Radiat. Trans- 
fer 

^ To whom all correspondence should be addressed. Phone: (+7) 812/428 42 63; 
Fax: (+7) 812/428 71 29; e-mail:nvv@aisbpu.spb.su 



Preprint submitted to Elsevier Preprint 



1 February 2008 



1 Introduction 



In many scientific and engineering applications prolate and oblate spheroids 
are appropriate models for real particles. The scattering of electromagnetic 
radiation by spheroids is usually calculated by the separation of variables 
method (0; D (SVM) and the T-matrix method (g) (TMM). Although both 
approaches are well-suited for numerical computations, it is often difficult to 
give reliable error estimates. For non-absorbing particles, there always exists 
an internal check as the absorption cross-sections (and efficiency factors) must 
be equal to zero. Obviously, for particles consisting of absorbing materials 
like silicates, graphite or metals, this test does not work. Despite the fact 
that rather elaborated methods for the calculation of the optical properties 
of highly absorbing non-spherical particles can be found in literature (see, for 
example, Michel et al. @), these methods have never been used for detailed 
numerical studies. 

In this paper, we consider the light scattering by homogeneous spheroids using 
the most popular approaches - SVM and TMM. The extinction and polariza- 
tion cross-sections for prolate and oblate spheroids of various aspect ratios a/b 
are calculated for the refractive indices m typical of soot, graphite, iron, sili- 
cate, and water ice and compared with the results for the particles of the same 
volume. A part of the results is presented in tabular form serving as bench- 
marks for forthcoming calculations. We also describe a combined approach to 
calculate the light scattering by axisymmetric particles. The method should 
allow to treat easily (very) elongated and flattened absorbing particles. 



2 Description of methods 

2.1 General definitions 

A spheroid (ellipsoid of revolution) is obtained by the rotation of an ellipse 
around its major axis (prolate spheroid) or its minor axis (oblate spheroid). 
The ratio of the major semiaxis a to the minor semiaxis b (i.e. the aspect ratio 
a/b) characterizes the particle shape which may vary from a nearly spherical 
one {a/b ^ 1) to a needle or a disk {a/b 3> 1). 

We assume an incident plane wave having the wavelength A. Let a denote the 
angle between the propagation direction and the rotation axis of the spheroid 
(0° < a < 90°). 

For axial propagation {a = 0°), there is no polarization of transmitted radi- 
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ation due to symmetry. If a ^ 0°, two cases of polarization of the incident 
radiation have to be considered: the electric vector E is parallel (TM mode) 
or perpendicular (TE mode) to the plane defined by the spheroid's rotation 
axis and the wave propagation vector. 



The size parameter is given by 



27rrv 

= (1) 



where ry is the radius of the sphere whose volume is equal to that of the 
spheroid. The radius ry is defined as 

r\ — a}? for prolate spheroids (2) 

and 

r\ = a% for oblate spheroids. (3) 

One usually calculates the efficiency factors Q = C /G which are the ratio of 
the corresponding cross-sections C to the geometrical cross-section G of the 
spheroid (the area of the particle's shadow) 

1 /2 

G{a) — nb (a^ sin^ a + b^ cos^ for prolate spheroids (4) 

and 

1 /2 

G{a) = na (^a^ cos^ a + 6^ sin^ aj for oblate spheroids. (5) 

In order to compare the optical properties of the particles of different shapes 
it is convenient to consider the ratios of the cross-sections for spheroids to the 
geometrical cross-sections of the equal volume spheres, G/nry. They can be 
found as 

G _[{a/by sin'' a + cos''a\^ 



for a prolate spheroid and 

G _ [(a/6)^cos^« + sin^c^]V^ 
(a/6)V3 

for an oblate spheroid. 



Q (7) 
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The albedo of a particle can be calculated from the extinction and scattering 
cross-sections 

^ sea 



c, 



ext 



The dichroic polarization efficiency is defined by the extinction cross-sections 
for TM and TE modes 



■ ■ 100%, (9) 



ext ~_ext 

vrM~T~7vrE 



^ext "T ^ext 



where the upper (lower) sign is related to prolate (oblate) spheroids. This 
ratio describes the efficiency of the polarization of light transmitted through 
the uniform slab consisting of non-rotating particles of the same orientation. 

The optical properties of spheroidal particles can be determined by various 
methods of light scattering theory. Most frequently, the separation of variables 
method obtained by Asano & Yamamoto (0) and Farafonov (Q ^) and the T- 
matrix method (^; P) are used. The calculations can be also performed with 
the discrete dipole approximation (DDA) which can be applied to arbitrarily 
inhomogeneous and irregular particles, but requires large computation time 
and is less accurate than the other methods (see discussion in the papers 

e I))- 



2.2 Asano & Yamamoto 's solution for spheroids (SVMl) 



Historically it was the first approach to the light scattering problem for spheroids 
with a complex refractive index. The method is based on the solution to the 
Helmholtz equation in the spheroidal coordinate system. 

In their pioneering paper, Asano & Yamamoto (|l]) used the Debye potentials 
to describe the electromagnetic fields, which is similar to the Mie solution 
for spheres. The scattering coefficients then are bound in the infinite systems 
of the linear algebraic equations and can be found by solving the truncated 
systems. 

In our comparison, we made the calculations with the numerical code of Peter 
Martin (see (|; |10D). 
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2.3 Farafonov's solution for spheroids (SVM2) 

The principal distinction of this solution (see (|^; |lT])) from the previous one 
is the special basis for the representation of the electromagnetic fields - a 
combination of the Debye and Hertz potentials (i.e. the potentials introduced 
to solve the light scattering problem for spheres and infinitely long cylinders, 
respectively). 

The approach has an incontestable advantage for strongly elongated or flat- 
tened particles. In this paper, the most recent version of our numerical code 
has been used. It involves different methods to calculate the spheroidal func- 
tions with the automatic choice of the most appropriate method among them. 

2.4 Solution for spheroids by the T-matrix method (TMMl and TMM2) 

Besides SVM, another very popular approach to solve the light scattering 
problem for spheroids is the T-matrix method. It is mainly applied to the 
axisymmetric particles: finite cylinders, spheroids, Chebyshev particles. 

The main idea of the method is the expansion of the incident, internal and 
scattered radiation fields in terms of the vector spherical harmonic functions. 
Because of the linearity of the Maxwell equations, the relations between the 
expansion coefficients of the fields are linear and given by two matrices. Thus, 
solution consists of calculations of the matrix elements being the surface in- 
tegrals and inversion of a matrix. The advantages of the approach are simple 
compact codes and their high speed. Another advantage is also the ability to 
treat more the complex case of elastic waves, the potential applicability to 
particles of complex shape, and the possible analytical averaging over the par- 
ticle orientations within the expansion coefficients which can greatly enlarge 
the field of the TMM applications. 

We used the T-matrix codes of Barber and Hill (|^) (TMMl) and Mishchenko 
(TMM2; see (^ and references therein). 

2.5 A comparison of methods 

We do not intend here to compare in detail the computational efficiency of 
different methods. The goal of this subsection is only to illustrate it on one 
example. A crude impression of the range of applicability of the methods may 
be obtained from Fig. 1 where the normalized cross-sections are plotted for 
oblate spheroids with m = 1.7 -|- 0.7i and a/b = 4. In all four cases double 
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Fig. 1. Normalized extinction cross-sections vs xy for oblate spheroids with the 
refractive index m = 1.7 + 0.7i, a/b = 4 and a = 0°. The solutions of Asano & 
Yamamoto (SVMl) and Farafonov (SVM2) and the T-matrix codes of Barber & 
Hih (TMMl) and Mishchenko (TMM2) were used in the calculations. Note that 
because of the special normalization (see Eq. (7)) the cross-sections have the limit 
2(a/6)2/3 PS 5.040 if 

precision was used. 

It should be noted that the advantage of the T-matrix codes appears mainly 
for particles more spherical than presented in the Fig. 1. On the other hand, 
for larger aspect ratios, the SVM2 code becomes unrivalled. 



2.6 A new method for axisymmetric particles 



Each of the methods (SVM and TMM) has its own problems. We suggest the 
following approach that could solve some of them. 

The incident, scattered and internal fields are divided in two parts: axisym- 
metric and non-axisymmetric ones 



E = El + E2, H = Hi + H2 



(10) 



It is possible to do so that the axisymmetric part does not depend on the 
azimuthal angle {(f) and the averaging of the non-axisymmetric part over all 
the azimuthal angle values gives zero (0). Then the scattering problem can be 
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solved for each of the parts separately. 

The scattering problem is formulated in the integral form like in the TMM 
(e.g., (^). The Abraham's potentials 

p = Ei^cosv2, q = Hi^cos(p (11) 



are introduced for the axisymmetric part. A superposition of the Debye and 
Hertz potentials is used for the non- axisymmetric part: 

E2 = V X {Ue, + Vr), H2 = x V x {Ue, + Vf) (12) 



for TE mode and 



ieki 



-V X V X (t/e, + yr), H2 = V x {Ue, + Vf) (13) 



for TM mode. Here e is the complex permittivity, n the magnetic permeability, 
ko the wavevector in vacuum, r the radius- vector, the unit vector, and i the 
complex unity. 

The scalar potentials are expanded in terms of the spherical wave functions. 
The coefficients of the expansions are determined from solution to the algebraic 
systems similar to those obtained in the TMM (e.g., (|^)). 

Thus, the new approach should combine the strong aspects of the basic meth- 
ods: the simple solution scheme typical of the TMM and the ability to treat 
particles whose shape can be very elongated or flattened as in the SVM2. The 



theoretical description of this approach is given in ([12D the development of a 
numerical code realizing it is in progress (P^). 



3 Numerical results 



Here we present some results illustrating the behaviour of the optical proper- 
ties of spheroidal particles for the case of transmitted radiation. Because this 
is the first consideration of optical properties of highly absorbing non-spherical 
particles beyond the Rayleigh limit, we try to demonstrate the features of such 
particles in extinction, scattering and polarization. 

Note that in the standard case of light scattering by spheres (Mie theory) 
only the refractive index and the size parameter of particles may be varied. 
For the homogeneous spheroids, one can also study the effects of the particle 
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shape (the aspect ratio a/h) and orientation (the angle a) as well as that of 
the polarization state of the incident radiation (TM and TE modes) for both 
elongated (prolate spheroids) and flattened (oblate spheroids) particles. 

We consider spheroids with aspect ratios a/h — 1.1,2,4, and 10 in a fixed 
orientation. The first a/h value illustrates the influence of small deviations 
from the spherical shape. All the calculations presented in this Section were 
performed using our new SVM2 code. 



3. 1 Eodinction 

Figures 2 and 3 show the normalized extinction cross-sections for non-absorbing 
and absorbing spheroids with different aspect ratios. The similar scale of the 
y— axes allows to see clearly the shape effects for the particles of equal volume. 

Note that large oscillations and the ripple-hke structure disappear for ab- 
sorbing spheroids - the curves become very smooth. These features depend 
mainly on the value of the imaginary part of the refractive index. The curves 
Cext /^rry (x) for a = 0° and a = 90° differ more and more when the ratio 
a/h grows (Fig. 3). For absorbing spheroids, the position of the first (main) 
maximum slightly shifts to larger values of xy with increasing a/h, while for 
non- absorbing particles, we have the opposite case. 

Now let us consider the particle of a fixed volume and begin to change its 
shape from a sphere to a needle or disk. The extinction by the particle is 
shown in Fig. 4 for three values of the refractive index. As follows from this 
Figure, in almost all cases the extinction cross-sections seem to reach their 
asymptotic values already at a/h ~ 10. It is interesting that the extinction 
by non-spherical transparent particles is usually smaller than that by spheres. 
The absorbing spheroids normally remove from the incident beam 1.5-2 
times more energy than spheres of the same volume. Though the behaviour 
of Cext is refractive index/size dependent, the conclusions we made are rather 
general for absorbing spheroids. From Figs. 3, 4 one can also see that even 
after averaging over particle's orientations we should get a stronger extinction 
for particles with larger aspect ratios a/h. 



3.2 Scattering and albedo 

If the imaginary part of the refractive index is not zero, we can study in- 
dependently the scattering and absorbing properties of particles. In the case 
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SPHEROIDS : prolate, m = 1.3 + 0.01 
a/b = 1.1 


a/b = 2 

moe^ 1 1 1 1 1 1 1 




- a/b = 4 


1 1 1 1 1 1 1 1 


QOOQO = 90". TM-mode 
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/\ 



SPHEROIDS : oblate, m = 1.3 + 0.01 
a/b = 1.1 




a/b = 3 




a/b = 4 




Fig. 2. Normalized extinction cross-sections vs for prolate and oblate spheroids 
with m = 1.3 + O.Oi 

presented as an example in Fig. 3, the behaviour of extinction is mainly de- 
termined by that of absorption. 

The integral scattering properties of particles are characterized by their albedo 
(see Eq. (8)). This quantity depends on the particle size and, in general, on 
the particle shape. As far as we know the shape effects have not yet been 
discussed in the literature. Figure 5 shows the size dependence of the albedo 
for spheroids with m = 1.7 + 0.7i. The calculations were made for the particles 
with four aspect ratios a/b — 1.1, 2, 4, 10 and for two orientations a — 0° and 
90° (we adopt that the incident radiation is non-polarized). 

It can be clearly seen that the albedo for non-spherical particles whose sizes 
are larger than a critical value {xy 2 — 3) does not deviate more than ~ 20 % 
from that of spheres. The largest deviations in Fig. 5 occur for the aspect ratio 
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Fig. 3. Normalized extinction cross-sections vs xy for prolate and oblate spheroids 
with m = 1.7 + 0.71 



a/b = 10 and a = 0° (prolate spheroids) and a = 90° (oblate spheroids). Note 
that for arbitrarily oriented particles (3D-alignment), the distinction between 
the albedo of spheres and spheroids should be smaller (see Fig. 11 in (^ where 
the albedo for particles with m = 1.53 + O.OOSz and a/6 = 2 is plotted). 

Our calculations made for particles with different absorption show that the 
distinction of the albedo for spheres and spheroidal particles remains rather 
small (within ~ 20 %) if the ratio of the imaginary part of the refractive 
index to its real part k/n ^ 0.2 — 0.3, which corresponds to k ^ 0.3 — 0.4 if 
n = 1.7. Thus, we can conclude that the albedo of spheroidal particles from 
soot, graphite, iron very slightly depends on the particle shape. 
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Fig. 4. Normalized extinction cross-sections (relative to that for equal vol- 
ume sphere) vs a/h for prolate and oblate spheroids with xy = 3 and 
m = 1.3 O.Oi, 1.7 + O.Oi, and 1.7 + 0.7i. 



3. 3 Polarization 



If a volume contains aligned non-spherical particles, the initially non-polarized 
incident radiation will be partially polarized after passing the volume. The 
polarization degree can be used to estimate the alignment degree of the parti- 
cles. The simplest and at the same time extreme case of particles' alignment is 
the perfect alignment of non- rotating particles (picked fence orientation). The 
maximum polarization usually occurs when the major axis of the particles is 
perpendicular to the direction of the incident radiation (a = 90°). 

The previous calculations of the polarization (see (§ |10|)) did not deal with 
large, highly absorbing particles. However, such particles exist in nature and 
can be aligned under terrestrial and extraterrestrial conditions. 

The behaviour of the polarization efficiency P/r (see Eq. (9)) for non-absorbing 
and absorbing spheroids is shown in Fig. 6 for the case when a maximum po- 
larization is expected (picked fence alignment, a = 90°). It is clearly seen 
from this Figure that the relatively large particles produce no polarization — 
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Fig. 5. Albedo of spheroidal particles with m = 1.7 + 0.7i normalized relative to 
albedo of spherical particles with the same refractive index. The solid lines connect 
the points for a = 0°, the dashed lines connect the points for a = 90°. 

independent of their shape. For absorbing particles, it occurs at smaller xy 
values than for non-absorbing particles. Note also that the polarization re- 
versal takes place for disk-like particles. This effect seems to depend on the 
imaginary part of the refractive index as well. 



3.4 Benchmarks 



The results of computations for spheroids were previously published in the 



tabular form only in the papers of Voshchinnikov & Farafonov (|Tl|; 113; |T5|) 
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Fig. 6. Polarization efficiency vs xy for prolate and oblate spheroids with 
m = 1.3 + O.Oi and 1.7 + 0.7i, a = 90°. 

(efficiency factors for spheroids with the refractive index m = 1.2), Kuik et 
al. (|16D and Hovenier et al. (0) (the scattering matrices for spheroids with 
m = 1.5+ki,k < 0.01). Extensive calculations of the elements of the scattering 
matrices for randomly oriented slightly absorbing spheroids with a/b <2 were 
performed by Mishchenko and others (see (^ and references therein). 

In Tables 1 - 4, we present the normalized extinction and scattering cross- 
sections for prolate and oblate spheroids with the refractive indices m = 2.5 + 
O.Oi, 1.7 + 0.7i, 2.5 + 1.5i and 3.0 + i.Oi. The results — all 7 (or in a few cases 
6) significant digits — were obtained with at least two of four codes used (see 
Sects. 2.2-2.4). 



Editor: Tables 1-4 are put here!!! 



The numbers printed by italics in Table 1 were calculated using the SVM2 code 
only. As the extinction and scattering cross-sections were the same (with 8 and 
more digits), the values presented in the Table are expected to be correct. 



4 Conclusions 



On the basis of the solutions to the light scattering problem by the separa- 
tion of variables method and the T-matrix method we have considered the 
optical properties of absorbing spheroidal particles of different aspect ratios. 
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Distinct approaches allowed to get the reliable results which may be used as 
benchmarks. 

The following conclusions can be made from the study: 

1. The extinction cross-sections for highly absorbing spheroids are typically 
1.5-2 times larger than those for spheres of the same volume. 

2. The albedo of non-spherical particles exhibits only a weak dependence on 
the particle shape and is determined mainly by the imaginary part of the 
refractive index. 

3. For particles larger than a certain minimum size, the spheroidal particles 
do not polarize the transmitted radiation independent of their shape. 
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Table 1: Normalized scattering cross-sections Csca./'^^v for spheroids with m = 2.5 + O.Oi 



Prolate spheroid Oblate spheroid 



a/b 


Xv 


a = 0° 


a = 90° (TM) 


a = 90° (TE) 


a = 0° 


a ^ 90° (TM) 


a = 90° (TE) 


2 





1 


8.132300E-05 


2.262901E-04 


8.162806E-05 


1.641502E-04 


5.756336E-05 


1.637629E-04 







3 


6.610480E-03 


2.018734E-02 


6.819180E-03 


1.433111E-02 


4.685330E-03 


1.404947E-02 




1 





1.299951E-00 


3.736834E-00 


9.854978E-01 


2.464760E-00 


6.647633E-01 


3.680223E-00 




3 





2.756951E-00 


2.423317E-00 


1.456577E-00 


4.484763E-00 


1.777398E-00 


3.846601E-00 


4 





1 


6.871662E-05 


4.258699E-04 


6.952952E-05 


2.612223E-04 


3.688507E-05 


2.600006E-04 







3 


5.105309E-03 


3.919277E-02 


5.615168E-03 


2.307084E-02 


2.874841E-03 


2.217663E-02 




1 





7.055768E-01 


5493155E-00 


5.42674OE-OI 


3.362340E-00 


2.595873E-01 


5.360019E-00 




3 





1.673046E-00 


5.7558UE-00 


4.730299E-00 


739 2485 E- 00 


2.580128 E- 00 


3.791884E-00 


10 





1 


6.058991E-05 


6. 745 045 E- 04 


6.315583E-05 


4.4235O6E-O4 


2.643817E-05 


4.383382E-04 







3 


3.U6015E-03 


5.950350E-02 


4.588377E-03 


3.852378E-02 


1.918230E-03 


3. 569 47 OE- 02 




1 





3. 0883 HE- 01 


6.153316 E- 00 


2.283901E-01 


3. 605985 E- 00 


1. 053865 E- 01 


5.126838E-00 




3 





1612520E-01 


9.781859E-00 


7081355 E- 00 


7 766600E-00 


7497780E-00 


4.7OI476E-OO 
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Tabic 2: Normalized extinction and scattering cross-sections Cext,sca/7rry for spheroids 
with m = 1.7 + 0.7i 



a = 0° a = 90° 



a/b 


Xv 




Cext/TI-r-y 




Csca/7rr^ 


































Prolate spheroid 
















2 


0, 


.1 


9, 


.260996E-02 


6 


.520100E- 


■05 


1.867292E-01 


1 


.323250E- 


■04 


9, 


.250492E-02 


6, 


.544660E-05 




0, 


.3 


3 


.GG1614E-01 


5 


.177635E- 


■03 


6.027228E-01 


1 


.099226E- 


■02 


2, 


.974687E-01 


5, 


.348186E-03 




1, 


.0 


1, 


.697352E-00 


3 


.846683E- 


01 


2.477888E-00 


7, 


.841288E- 


■01 


1, 


.556310E-00 


4, 


.492447E-01 




3 


.0 


2. 


.481402E-00 


1 


.039641E- 


■00 


3.195685E-00 


1 


.600695E- 


■00 


3 


.021085E-00 


1. 


.360101E-00 




10 


.0 


1, 


.8G433GE-00 


9 


.737252E- 


■01 


2.874307E-00 


1 


.607439E- 


■00 


2, 


.908106E-00 


1, 


.453882E-00 


4 


0, 


.1 


8, 


.155853E-02 


5, 


.681280E- 


■05 


2.550670E-01 


1 


.807060E- 


■04 


8, 


.14151 E-02 


5, 


.748690E-05 




0, 


.3 


2 


.62113 E-01 


4, 


.141768E- 


■03 


8.11370 E-01 


1 


.476105E- 


■02 


2, 


.585912E-01 


4, 


.568657E-03 




1, 


.0 


1, 


,50146 E-00 


2 


.669788E- 


■01 


2.713089E-00 


8 


.506337E- 


■01 


1, 


.24346 E-00 


3 


.192207E-01 




3, 


.0 


2, 


.580295E-00 


9 


.865859E- 


■01 


3.848184E-00 


1 


.966160E- 


■00 


3 


.40958 E-00 


1, 


.515164E-00 


10 


0, 


,1 


7, 


.60280 E-02 


5, 


.084200E- 


■05 


3.01197 E-01 


2, 


.122690E- 


■04 


7, 


.58943 E-02 


5, 


.299880E-05 




0, 


.3 


2, 


.39882 E-01 


2 


.843924E- 


■03 


9.378379E-01 


1 


.631801E- 


■02 


2. 


.368506E-01 


3 


.807946E-03 


















Oblate spheroid 
















2 





,1 


1, 


.533476E-01 


1 


.086960E- 


■04 


6.964599E-02 


4 


.913110E- 


■05 


1, 


.534705E-01 


1, 


.084390E-04 




0, 


.3 


4, 


.941085E-01 


9 


.029149E- 


03 


2.246184E-01 


3 


.938983E- 


03 


4, 


.971749E-01 


8, 


.843852E-03 




1, 


.0 


2, 


.198062E-00 


7, 


.105949E- 


■01 


1.291186E-00 


3 


.151001E- 


■01 


2, 


.326983E-00 


6, 


.221156E-01 




3, 


.0 


3 


.378955E-00 


1 


.672422E- 


■00 


2.603296E-00 


1 


.040019E- 


■00 


2, 


.768031E-00 


1, 


.330862E-00 


4 


0, 


.1 


2, 


.021662E-01 


1 


.432990E- 


■04 


4.737106E-02 


3 


.331120E- 


■05 


2, 


.023435E-01 


1, 


.426270E-04 




0, 


.3 


6, 


.442415E-01 


1 


.176912E- 


■02 


1.50949 E-01 


2 


.574738E- 


■03 


6, 


.483922E-01 


1, 


.129677E-02 




1, 


.0 


2, 


.39421 lE-00 


7, 


.756880E- 


■01 


8.254487E-01 


1 


.707659E- 


■01 


2. 


.65821 lE-00 


6. 


.204606E-01 


10 


0, 


.1 


2, 


.59201 E-01 


1 


.833420E- 


■04 


3.51934 E-02 


2 


.458170E- 


■05 


2 


.59355 E-01 


1, 


.816740E-04 




0, 


.3 


8, 


.101747E-01 


1, 


.452893E- 


■02 


1.102168E-01 


1 


.774836E- 


■03 


8, 


.13450 E-01 


1, 


.343542E-02 
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Table 3: Normalized extinction and scattering cross-sections Cext,sca/7rry for spheroids 
with m — 2.5 + 



a = 0° a = 90° 

a/b XV C^^/nr'y C^^J^ ~C™/^ C^J^ C^J^ C^J^ 



Prolate spheroid 



2 


0, 


.1 


7, 


.015863E-02 


1 


.324916E- 


■04 


2.555773E-01 


4, 


.901860E-04 


6, 


.978007E-02 


1 


.329890E- 


■04 




0, 


,3 


2, 


.606615E-01 


1 


.085718E- 


■02 


9.351659E-01 


4, 


.418671E-02 


2, 


.499507E-01 


1 


.120452E- 


■02 




1, 


,0 


2, 


,545641E-00 


8 


.269798E- 


-01 


4.506818E-00 


2, 


.290117E-00 


2, 


.378954E-00 


1 


.018659E- 


-00 




3, 


,0 


2, 


.543020E-00 


1 


.196810E- 


-00 


3.377612E-00 


2, 


.030276E-00 


3 


.018778E-00 


1, 


.561478E- 


-00 


4 


0, 


.1 


5, 


.77412 E-02 


1 


.079560E- 


-04 


5.61456 E-01 


1, 


.078514E-03 


5, 


.73428 E-02 


1 


.092344E- 


-04 




0, 


,3 


2, 


.09854 E-01 


8 


.054995E- 


-03 


1.994433E-00 


9, 


.529425E-02 


1, 


.99005 E-01 


8 


.869991E- 


-03 




1. 


.0 


2, 


.23898 E-00 


5 


.318469E- 


■01 


4.844781E-00 


2, 


.443652E-00 


1, 


.709826E-00 


6 


.654597E- 


■01 




3. 


.0 


2, 


.564187E-00 


9 


.810827E- 


■01 


4.17736 E-00 


2. 


.576599E-00 


3 


.444722E-00 


1 


.733285E- 


■00 


10 


0. 


.1 


5, 


.18757 E-02 


9 


.344730E- 


■05 


8.83893 E-01 


1, 


.689619E-03 


5, 


.15928 E-02 


9 


.740630E- 


■05 


















Oblate spheroid 
















2 


0, 


.1 


1, 


.683208E-01 


3 


.228280E- 


-04 


4.590543E-02 


8, 


.685830E-05 


1, 


.690525E-01 


3 


.220660E- 


-04 




0, 


,3 


6, 


.046926E-01 


2, 


.854920E- 


-02 


1.682741E-01 


7, 


.112129E-03 


6, 


.24682 E-01 


2, 


.798368E- 


-02 




1, 


,0 


4, 


.007697E-00 


2 


.005398E- 


■00 


1.852325E-00 


6, 


.182319E-01 


3 


.928646E-00 


1 


.726163E- 


■00 




3, 


,0 


3, 


.644115E-00 


2 


.139153E- 


-00 


2.436512E-00 


1, 


.057701E-00 


2, 


.916411E-00 


1 


.654727E- 


-00 


4 


0, 


.1 


3, 


.07216 E-01 


5 


.903690E- 


■04 


2.714797E-02 


5, 


.130230E-05 


3, 


.085216E-01 


5, 


.876060E- 


■04 




0, 


,3 


1, 


.089947E-00 


5, 


.221910E- 


-02 


9.612426E-02 


4, 


.009507E-03 


1 


.123093E-00 


5, 


.017062E- 


-02 




1. 


.0 


4, 


,.1.10.19<SE-()() 


2 


.;-!l;-!7(i7E- 


-00 


i.04;-')0;-') E-00 


2 


,S7(i447E-01 


4 


,()S404 E-00 




.<S474<S<SE- 


-00 
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Table 4: Normalized extinction and scattering cross-sections Cext,sca./'^fv spheroids 
with m = 3.0 + 4.0i 



a 



0° 



a 



90° 



a/b XV Cext/7rr^ 



Csca/7rr^ 



Csca /■"■''V' Qxt/ 



Csca/ 



0.1 3.459049E-02 

0.3 1.944302E-01 

1.0 1.894857E-00 

3.0 2.388466E-00 

0.1 1.017897E-01 

0.3 4.583023E-01 

1.0 5.236749E-00 

3.0 3.647681E-00 



1.836889E-04 
1.537309E-02 

9.206887E-01 
1.478325E-00 

5.864570E-04 

5.527781E-02 
3.731015E-00 
2.814034E-00 



Prolate spheroid 

1.929191E-01 1. 

9.098178E-01 1. 

5.495079E-00 3. 

3.453263E-00 2. 

Oblate spheroid 
2.106561E-02 1. 
1.197910E-01 9. 
1.404715E-00 6. 
2.009705E-00 1. 



115412E-03 
129382E-02 
970900E-00 
702600E-00 

116025E-04 

258967E-03 
821323E-01 
173967E-00 



3.339258E-02 
1.661101E-01 

2.283547E-00 
2.762240E-00 

1.046258E-01 

5.167953E-01 
4.215200E-00 
3.040199E-00 



1.843746E-04 
1.585916E-02 

1.438695E-00 
1.941164E-00 

5.850970E-04 

5.426552E-02 
2.800509E-00 
2.239816E-00 
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